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Abstract 

We study a previously introduced bi-local gauge invariant reformulation of 
two dimensional QCD, called 2d HadronDynamics. The baryon arises as a 
topological soliton in HadronDynamics. We derive an interacting parton model 
from the soliton model, thus reconciling these two seemingly different points of 
view. The valence quark model is obtained as a variational approximation to 
HadronDynamics. A succession of better approximations to the soliton picture 
are obtained. The next simplest case corresponds to a system of interacting 
valence, 'sea' and anti-quarks. We also obtain this 'embellished' parton model 
directly from the valence quark system through a unitary transformation. Using 
the solitonic point of view, we estimate the quark and anti-quark distributions of 
2d QCD. Possible applications to Deep Inelastic Structure Functions are pointed 
out. 
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1 Introduction 



There are currently two distinct points of view on what a baryon is. One may 
be traced back to the quark model. In the other point of view, baryons arise as 
solitons of low energy effective local field theories of mesons, an idea that may 
be traced back to the Skyrme Model Q^. In previous work by one of us, it was 
shown that in two dimensions, there is an exact description of QCD as a bi-local 
theory of mesons This description, called 2 dimensional HadronDynamics, 
is not a low energy effective theory, but is equivalent to 2d QCD for all energies 
and numbers of colors. The exact description of the baryon is as a topological 
soliton of this bi-local theory. At the other extreme, we studied a 2 dimensional 
interacting quark model for the structure of the baryon 3 . The question then is 
whether the quark model picture can be derived from the exact bilocal solitonic 
picture in two dimensions. Here we derive the interacting valence quark model 
as a variational approximation to this bi-local soliton theory. Moreover, we find 
a succession of increasingly accurate variational approximations to the soliton 
model. The next simplest case turns out be a system of interacting valence, 
'sea' and anti-quarks. We show that this 'embellished' parton model could 
also have been obtained directly via a unitary transformation applied to the 
valence quark model, recovering the 'Bogoliubov' transformation introduced in 
Thus we have a reconciliation of the exact bi-local soliton model with the 
simpler relativistic parton picture of the baryon, in 2 dimensions. 

The main advantage of this new point of view is that the semi-classical 
approximation of HadronDynamics corresponds to the large Nc limit of 2d QCD, 
and so is capable of describing non-perturbative phenomena such as the struture 
of hadrons. Thus this reconciliation between the soliton and parton pictures is 
more than just a mathematical correspondence. We illustrate its usefulness 
by calculating approximately, the quark and anti-quark distributions of the 
baryon in 2d QCD. In fact, due to the correspondence we establish here, the 
approximate solutions of the quark models presented in jSj are actually also 
approximate solutions of 2d QCD. Thus, in this paper, we will focus more on 
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the passage from the bi-local sohton theory to quark models, rather than their 
actual solution. As mentioned in 3 , our results agree well with the direct 
numerical solutions of Hornbostel et. al. 4 . To summarize, we find that the 
valence quark approximation is accurate not only in the non-relativistic limit, 
but also in the ultra-relativistic chiral limit. In particular, we find that in the 
chiral limit, the valence quark approximation is exact, for Nc — s- oo. 

As an aside, we speculate on the possible phenomenological implications of 
the above model. In Deep Inelastic Scattering [S], the transverse momenta of 
the partons is small compared to their longitudinal momenta. Moreover, the 
observables of interest, the baryon structure functions, depend only on the pro- 
ton momentum P and photon momentum q, which lie in a two-dimensional time 
like hypersurface, spanned by time and the beam direction. Thus, there must 
be an effective 2 dimensional theory that describes the structure of the proton 
as measured in Deep Inelastic Scattering. Two dimensional HadronDynamics 
has the correct symmetries to be a candidate for an approximate description of 
the 'relevant' interactions of the quarks, in such an effective action. This model 
can be thought of as a representative example of such an approximate effective 
action. The variational and many-body techniques developed here should be 
useful in understanding any such 2 dimensional effective action. 

2 Two Dimensional Quantum Hadron Dynam- 
ics 

Let us begin with a summary of Two dimensional Quantum HadronDynamics 
12]. 2d QCD is quantized in the null gauge = in null coordinates. The 
elimination of longitudinal gluons leads to a linear potential between quark 
fields a, which satisfy canonical anti-commutation relations. The resulting 
hamiltonian is: 

^ = /dxat"i[p + If]aa^ - i / - y| : at«(a;)a„,(x) :: a^^Hy)at^{y) : dxdy. 
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'g' is a coupling constant with the dimensions of mass, i,j are color indices and 
a, b are flavour indices. Define the bilocal gauge invariant variable M^{x,y) = 
: au{x)a^°''^{y) :. The points x,y are null separated. The operators M^{x,y) 
form a complete set of observables in the color singlet sector of two dimen- 
sional QCD. They provide a (projective) unitary irreducible representation of 
the infinite dimensional unitary Lie Algebra: 

[mI{p, q), M,{r, s)] = ^^{5l2i,5{q ~~ rWd sgn [p - s) + M.ip, s)]) 
(1) _-l(<5-27r^(s - pM sgn {r - q) + M,{r,q)]). 

Note that the commutators are of order -j^ so that the large limit is 
a sort of classical limit: plays the role that % does in an ordinary field 
theory. In this classical limit the above commutators are replaced by the Poisson 
Brackets of a set of classical dynamical variables M^{x,y). It was shown that 
the phase space of this system is an orbit of the unitary group, an infinite 
dimensional Grassmannian Gri . It is the set of all hermitean operators M with 
integral kernel M^{x,y) satisfying the quadratic constraint [e + M]^ = 1, with 
J \M{x,y)\'^dxdy < oo. Gri is a curved manifold with connected components 
labelled by an integer. The quadratic constraint is just the Pauli principle for 
quarks: the density matrix (p = i(l — M — e)) must be a projection operator. 
Here e is the Hilbert transform operator, diagonal in momentum space, with 
e{p,p) = sgn (p). 

So in the large Nc limit, our problem reduces to solving the equations of 
motion obtained from the hamiltonian and Poisson Brackets: 

(^^ = ~i/^ + 7"^^«^^'^^S + y / Ki^,y)M'a{y,x)\x-y\dxdy 

l{M,«(a;, y), M^{z, u)} = 5t5{y - z)[e2(x, u) + M^{x, u)\ 
(3) -d':,d{x-u)[el{z,y) + M-^{z,y)]. 
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The parameter /i^ is related to the current quark masses nia by a finite renor- 
maUzation: /i^ = m^ — ^ and = g^N^.- Though the hamiltonian is quadratic, 
this is a non-linear interacting theory since the phase space is a curved manifold 
due to the constraints on M[x, y). The linearization of the equation of motion 
around the vacuum A'l^ — describes an infinite number of free mesons, and 't 
Hooft's integral equation jH] for the meson masses was recovered. 

What kind of solution to this theory represents the baryon? The quantity 
was shown to be an integer, a topological invariant and 
hence conserved under time evolution. We sec that B is in fact baryon number. 
Thus, for p > 0, —\M^{p^p) and —^M^{—p, —p) represent the quark and anti- 
quark probability densities in the baryon. Thus the baryon is a topological 
soliton in this picture. It corresponds to a static solution of the equations of 
motion (minimum of energy subject to constraints), that has baryon number 
one. A Lorentz invariant formulation is to minimize the (mass)^ of the baryon: 

-i / Mip,p)^^ + ij dxdy\M{x, 2/)pi|x - y\ 

3 Separable or Rank One Ansatz and Valence 
Quark Model 

We have developed a method pi to find the minimum of energy on the the phase 
space: a variant of the steepest descent method that takes into account the 
non-linear constraint. Here, we describe another method based on variational 
approximations, which brings out the connection to the quark model. The 
main difficulty in minimizing the energy is the non-linear constraint satisfied by 
M{x,y). We find a succession of variational ansatzes for M{x,y) that replace 
this constraint with simpler ones. These define an ascending family of sub- 
manifolds, which form a dense subset of the phase space. Minimizing the energy 
on these sub-manifolds will give us successively better approximations. These 
variational ansatzes turn out to correspond to interacting quark models. 
To start with, consider an ansatz of the separable form 
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(4) MS{p,q) = -2r{p)rbiQ)- 

This satisfies the constraint (e + M)^ = 1 if i/! is of norm one and of positive 
momentum. The Poisson Brackets of the M^{x,y) imply the relations 



(5) {%{p),M<l)} = = {r*{p)A'*{q)}, 

(6) {Mp),r\q)} = -i27rSl6{p-q). 

The 'ipa by themselves define a classical dynamical system with hamiltonian 

We can quantize this 'mini' theory by looking for operators satisfying canonical 
commutation relations. Let us denote the parameter that measures the quan- 
tum correction, analogous to h, by The constraint on the norm can be 
implemented by restricting attention to those states \ V > satisfying 

(8) I ^ (p)^„(p)^|y>=i. 

A representation for our commutation relations is provided by bosonic creation- 
annihilation operators: 



(9) ikip), k{p')] = = [b\p), b\p')] , [k{p), b\p')] = 27r5(p - q)Sl 
with 

(10) ^a(x) = -^kix) , ^ta(a:) = -^i^^x)- 

Then the constraint becomes the condition that we restrict to states con- 
taining Nc particles: b'^°'{p)ba{p)^ = N^- Therefore, must be a positive 
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integer! Thus we are dealing with a system of Nc bosons interacting through a 
Hnear potential. 

What are these bosons? They are the valence quarks of the parton model. 
They appear like bosons in the momentum and spin-flavour quantum numbers 
since their wave function is totally anti-symmetric in color. Nc is interpreted as 
the number of colors. In the mean field approximation their wave function 
is e'i'"*™<^ V'(Pi) ■ ■ ■ '4'{PNa)i which corresponds to 

(11) \V>=a'l--af^\0>. 

If /3^(p, q) = -^a^°'''{p)abi{q) is the quark density operator, then the expectation 
value of p{p,q) in the mean field state state \V > is equal to the 'classical' 
density matrix pi{p,q) = ^{S{p,q) — Mi{p,q) — e{p,q)). Thus the classical (or 
large Nc) limit we have been discussing is just the mean field approximation 
to this many-body problem, an idea that goes back to Witten [7]. The semi- 
classical approximation will give us the leading corrections in the case of finite 

Since the quark null momenta are positive, their sum must equal the total 
baryon momentum P. In particular, the parton momenta cannot exceed P. 
This ensures that the quark distributions vanish beyond p = P. However, the 
total baryon momentum is extensive P ^ Nc- So in the limit as Nc oo, 
< p < oo. The valence quark wavefunction is determined by minimizing the 
{massY' subject to the normalization and momentum sum rule conditions: 



/o'^iv^wPH^i, Nc!^p\m?^^p- 

Here ^{x) — J^^ rp{p)e^P^ Since we have ignored them, we will get the spin 
and flavor averaged wavefunction. The probability density of valence quarks is 

2 



V{xb) 



p_ 

2-K 



^P{xbP) 



, where xg ^ ^. 

The ground state of this many body problem was found in [3 , in the guise of a 
valence quark model. Let us just summarize the result. In the limit Nc — s- oo and 
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m = 0, the absolute minimum of the variational principle is ip{p) — y^e^^p. 
Here P is the mean baryon momentum per color, It turns out that A4 = 
for this configuration, so that it is not just a minimum on the separable 
submanifold, but on the entire phase space, in the chiral and large Nc limits. 
Thus we find that the ground state baryon is massless in this limit, in agreement 
with Hornbostel et. al. 0]. 

A variational approximation to the ground state, after including the leading 
efltects of finite Nc was also given in In the chiral limit, the valence quark 
probability distribution is V{xb) = {N^ — 1)[1 — xb]^"^"^. This variational 
approximation agrees well with our numerical solution 3 and is identical to the 
numerical solution of Hornbostel et. al. (See ref. 4 Eqn. 22). 



4 Rank Three Ansatz: Valence, Sea and Anti 
Quarks 

We can get a better approximation to the exact soliton model, by considering a 
larger submanifold of the phase space, compared to the separable ansatz, which 
corresponded to the valence quark approximation. 

The departure from the valence quark picture is determined by the dimen- 

2 

sionless ratio a measure of chiral symmetry breaking. Thus we should 
expect the anti-quark content to be small for small current quark masses. The 
leading effect of finite Nc is to constrain the range of momenta of the partons, 
as we have seen. 

The mathematical advantage of the separable ansatz is that it 'solves' the 
nonlinear constraint on M: more precisely, it replaces it with the condition that 
tjj is of norm one. In the same spirit, consider the configuration 

r 

(12) Mr^ ^feV'a®^^'- 

a,b=l 

Here we choose tpa to be a set of r orthonormal eigenvectors of the operator 
e; i.e., e^pa = — il- This implies that the operator Mr is of rank r: 
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the special case of rank one is just the separable ansatz above. This ansatz will 
satisfy the constraint on M if the r x r matrix ^ is hermitean and satisfies the 
constraint 



(13) eaeb + [ea + ec]C = 0, 

a 'mini' version of the constraint on M. Moreover, the baryon number is -B = 
— |tr M = —^tr^. In the special case of rank one, we have simply ^ = —2. 
By choosing a large enough value of r this ansatz can produce as general a 
configuration in the phase space as needed. The simplest configuration of baryon 
number one that departs from the separable ansatz is of rank three. We will 
find that for small current quark masses, even this departure is very small, so 
we do not need to consider configurations of higher rank. 

By a choice of basis among the ipa, we can always bring a rank three config- 
uration of baryon number one to the form 



Ms = -2ti) (g) + 2C! [V'- (g) Vl - V'+ «) 
(14) +2C-C+[V'- V'l +-0+ ® -01-], 

where tp-j^j^^ are three vectors in L'^{R) satisfying eip^ = —V'-, = 
tp, eV+ = V'+.IIV'-IP = IIV'IP = II^+IP = 1, < '0,1/'+ >= 0. The conditions 
< ip-Tip >=< >= arc then automatic. The parameter < C- ^ 1 

measures the deviation from the rank one ansatz and liencc!, the anti-quark 



content of the baryon. (^_|_ = y 1 — For example, baryon number is given by 

(15) B = J^^ {mp)\'+c- [\Mp)f - i^-(-f)i']} f • 

V', '>p+ vanish for p < and describe valence and 'sea' quarks. Their orthogonal- 
ity is a consequence of the Pauli principle, ip- is the anti-quark wavefunction. 

2 

From our previous result we expect C- to vanish as ^ — > 0. 

This rank 3 ansatz can also be understood as arising from a unitary trans- 
formation applied to the valence quark ansatz. The phase space of HadronDy- 
namics carries a transitive action of the infinite dimensional restricted unitary 
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group PI- Thus the configuration M3 can be obtained from Mi by a unitary 
transformation 

(16) C/+(e-2V'(^ V^)f/ = e + e6V'a<»V'^^ 

Since both Mi and M3 have the same baryon number, U hes in the connected 
component of the identity and is of the form U = e''^ for A hermitean. From the 
above expressions for Mi , M3, we see that U is the identity except on the span of 
ip+ and ip-. e — — CT3, in this sub-space. Thus e~''^(— 0-3)6*"^ = sai + (r — 1)0-3, 
where r = 2Ci and s — 2(-C,+ and ai are the Pauh matrices. Therefore, on this 
subspace, A is a 2 * 2 traceless hermitean matrix cr.w. w is the vector in R'^ 
about which (0, 0,-1) must be rotated by an angle 2|w| to reach (s, 0, r — 1)). 
Thus A = i arcsin(C-)('0- -0+ - V'+ ® V'- ) • 

We can use the infinite dimensional analogue of the Pliicker embedding 011] 
of the Grassmannian in the Fermionic Fock space to reexpress this unitary trans- 
formation on the phase space of HadronDynamics, as a Bogoliubov transforma- 
tion on the second quantized states. The operator that corresponds to U and 
acts on the fermionic fock space is 

the sum over colors produces a singlet. The angle 9 of S can be identified 
as arcsin(C_). The second quantized state after the Bogoliubov transforma- 
tion is thus |FS'A >= e^'^^\V >. Here > is the valence quark state, and 
|FS'A > stands for a state containing valence, sea and anti-quarks. The condi- 
tion < VSA\p{p, q)\VSA >— p3{p, q) is then automatic, since the corresponding 
condition was satisfied in the rank 1 case and we have performed the same uni- 
tary transformation on both sides. 

Thus, we have derived the 'embellished' quark model, which contains va- 
lence, sea and anti-quarks, as a variational approximation to the bilocal soliton 
theory. The wave functions ij}, ip± and the probability of finding an anti-quark 
in a baryon . ~ a , were estimated in P]. To summarize, we found that the 
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probability of finding an anti-quark in the baryon is ~ 0.1% and a ^ .035 for 
Sj- ~ 10"'^. Moreover, the anti-quarks carry less than 0.1% of the baryon mo- 
mentum. Hornbostel et. al. 4 also find a similar suppression of the anti-quark 
content in the chiral limit. 

5 Conclusion 

Thus, while neither the quark model, nor the solitons of low energy effective local 
field theories provides a complete description of baryons in two dimensions, a bi- 
local quantum field theory provides an exact description. We have shown how 
the quark model arises as a variational approximation to this bi-local soliton 
theory. Moreover, we have used this bi-local theory to estimate the quark and 
anti-quark distributions in two dimensions and find good agreement with direct 
numerical approaches. It is interesting to know what the analogous non-local 
theory is in four dimensions. 
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